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ABSTRACT 

Galactic outflows of low ionization, cool (^ 10"* K) gas are ubiquitous in local starburst galaxies, and in the 
majority of galaxies at high redshift. How these cool outflows arise is still in question. Hot gas from supernovae 
has long been suspected as the primary driver, but this mechanism suffers from its tendency to destroy the cool 
gas as the latter is accelerated. We propose a modification of the supernova scenario that overcomes this 
difficulty. 

Star formation is observed to take place in clusters; in a given galaxy, the bulk of the star formation is found in 
the ^ 20 most massive clusters. We show that, for L* galaxies, the radiation pressure from clusters with Md > 
IQ^Mq is able to expel the surrounding gas at velocities in excess of the circular velocity of the disk galaxy. 
This cool gas can travel above the galactic disk in less than 2 Myr, well before any supernovae erupt in the 
driving cluster. Once above the disk, the cool outflowing gas is exposed to radiation, and supernovae induced 
hot gas outflows, from other clusters in the disk, which in combination drive it to distances of several tens to 
hundreds of kiloparsecs. Because the radiatively driven clouds grow in size as they travel, and because the hot 
gas is more dilute at large distance, the clouds are less subject to destruction if they do eventually encounter 
hot gas. Therefore, unlike wind driven clouds, radiatively driven clouds can survive to distances ^ 50kpc, and 
potentially give rise to the metal absorbers seen in quasar spectra. We identify these cluster-driven winds with 
large-scale galactic outflows. In addition, we note that the maximum cluster mass in a galaxy is an increasing 
function of the galaxy's gas surface density. As a result, only starburst galaxies, where massive clusters reside, 
are able to drive winds cold outflows on galactic scales via this mechanism. We find that the critical star 



formation rates above which large scale cool outflows will be launched to be S"' 
is in good agreement with observations. 

Subject headings: absorbers: Mg II - star formation rate - star cluster - outflows 
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1. INTRODUCTION 

Observational studies of both nearby and high redshift star 
forming galaxies have established that cold gas emerges from 
such galaxies at velocities ranging from a few tens to several 
hundr ed kilometers per secon d (Heckman, Armus, & Mile^ 
1990; Steideletal. 1996; Franxetal. 1997; Heckman etalJ 
2000; Pettini et al 2000 2001 ; Shaple v et al. 2003; Martin 
Rupke et al.ll2005HTrem onti et al.l l2007t" TWeiner et al.J 



2005 



2009t 



Menard et al.ll2009h . Hot outflowing gas from super- 



novae (SN) has been considered to be the primary mecha - 
nism for driving galactic outflows (IChevaUer & Clegil ( 119851) ; 
[Heckman (2002) and references therein). It is generally as- 
sumed that cold gas and dust are entrained in the hot gas. 

As numerous authors have pointed out, in such a sce- 
nario the cold gas i s subject to Kelvin- Helmholtz instabilities 
dKlein et al.lll994l; ICooper et al.l 120081) and/or by conductive 
evaporation dMarcolini et al.ll2005l) and ends up surviving for 
less than 1 Myr (a few cloud crushing times) before being de- 
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stroyed. The cloud crushing time follows fee = (a/pa)'''^ a/vb, 
where pd Ph is the density ratio between the cold gas and the 
hot gas, a is the radius of the cloud before it is exposed to 
the hot gas, and vi, w 1000km s~' is the shock velocity in the 
hot gas. This implies that the cold gas can travel less than a 
few hundred parsecs before being mixed into the hot gas and 
hence effectively destroyed. However, observational evidence 
shows that cold gas survives at least out to ten kpc. An ex- 
planation for the presence of this cold gas is required. In this 
paper we show how taking into account the effect of radiation 
pressure on dust grains can affect the fate of the gas. 

A second point regarding the driving of outflows is that 
most stars form in massive clusters, both in quiescent spirals 
that lack strong outflows, like the Milky Way, and in starburst 
galaxies such as M82. We argue that outflows emerge from 
the most massive clusters in a galaxy, as opposed to arising 
from the collective effects of smaller clusters. The radiation 
pressure from small clusters will not punch holes in the gas 
disk, and the subsequent supernovae will either leak into the 
interstellar medium (ISM) in low star formation rate galaxies, 
or will radiate their energy away after their growth has been 
halted by the high pressure of the ISM in high star formation 
rate galaxies (such as Arp 220). In contrast, the radiation pres- 
sure from massive clusters will blow cold gas out of the disk, 
paving the way for the hot gas from SN to escape. 

In this paper we present a model including both radiation 
pressure from massive clusters and ram pressure from super- 



novae. As recently pointed out by iNath & SiUd (l2009h . we 
find that before any cluster stars explode as SNe, radiation 
pressure, acting alone, launches the cold gas surrounding a 
star cluster to heights of several hundred parsecs above the 
disk. As the cold cloud travels above the disk, it expands. 
This expansion, as well as a lower hot gas density at larger 
radii, significantly increase the cloud crushing time once it 
does encounter hot gas, allowing it to survive long enough to 
reach distances exceeding one hundred kpc. Once above the 
disk, a cold gas cloud experiences a ram pressure force which 
is com parable to the radiation pressure force (jMurrav et al] 
lIOOl) . 

While it is clear that both driving mechanisms are signif- 
icant, the action of the radiation pressure is crucial for the 
survival of cold gas out to large radii. We also illustrate the 
importance of this force by showing that it is capable of driv- 
ing outflows even in the absence of hot gas. In ^we describe 
how the masses of the largest star clusters in a galaxy are de- 
termined by the interplay between self-gravity and radiation 
pressure in giant molecular clouds. In ^we describe how 
winds are launched from massive clusters, and how this leads 
to a critical star formation rate to produce galactic winds. We 
follow the evolution of these winds outside the star forming 
disk, to distance of order the galactic virial radius, in ^ We 
describe the cloud properties, cover factor, and mass loss rates 
of the wind in ^ We discuss our results in ^and offer con- 
clusions in the final section. 

2. PROPERTIES OF STAR CLUSTERS 

We show below that the cluster mass and radius determine 
the dynamics of the interstellar medium (ISM) around the 
cluster, out to distances comparable to the disk scale height. 
We describe how to estimate the cluster mass as a function of 
the surface density of the parent giant molecular cloud. We 
then describe the (observationally determined) cluster radius 
as a function of cluster mass. Finally, we describe the cluster 
mass distribution function, and conclude that most of the lu- 
minosity in a starburst galaxy is typically produced by a dozen 
or so massive star clusters. 



2.1. Cluster masses 

Star clusters form in dense cores inside GMCs. As shown 
by iMurrav et al.l ( 120101) . radiation from the young stars will 
transfer momentum to the surrounding gas through dust 
grains. The energy from the absorbed photons is then reemit- 
ted isotropically in the infrared; in most galaxies these IR pho- 
tons simply escape the galaxy. Q The bolometric luminosity 
of a young cluster is carried predominantly by ultraviolet ra- 
diation. The dust opacity in the ultraviolet is typically n ^ 
lOOOcm-g"', so these ultraviolet photons see an absorption 
optical depth of order unity at a column of Nh ~ 10~' cm"-. 
Since the gas in the vicinity of a massive proto-cluster has a 
density in excess of lO'^cm"'', the corresponding length scale 
is about 10'^ cm, i.e. 0.03 pc. As a result, tuv > 1. In the 
single-scattering regime (where t^/r < 1) essentially all pho- 
tons emitted by the cluster are absorbed once and then leave 

^ If the optical depth in the far-infrared, tfir, is larger than unity, which 
is the case in ultraluminous infrared galaxies (ULIRGs) and submillimeter 
galaxies, photons will experience multiple interactions with dust grains. This 
can enhance the effectiveness of radiation pressure by a factor of Tfm. 



the system as infrared photons. It follows that the radiation 
pressure force is 

^md^-. (1) 



The accretion of gas and the formation of stars occur un- 
til the radiation pressure becomes comparable to the gravi- 
tational force acting on the infalling material: Frad ~ fgrar- 
At this stage, the infall slows down significantly and likely 
takes place t h rough a disk. By equating these two forces, 
iMurrav et al.l (120101) showed that the fraction of gas ec = 
Md/Mc in a GMC that is converted into stars is proportional 
to the surface density Sc of the parent GMC, for GMCs that 
are optically thin to far-infrared radiation. 



£0 = 



ttGc , 

2{L/mS 



(2) 



where L/M^, is the light to mass ratio of the cluster This 
relation shows that the characteristic (or maximum) mass of 
the star clusters in a galaxy will scale with the surface density 
of the galactic disk S^ and hence the star formation rate of 
the hos t ga laxy. We will make use of the above relation in 
section [3T2I 



The star formation process is believed to be rather ineffi- 
cient. Typical star formation rates in massive (10^- IQ^Mq) 
proto-clusters are ^ 1 Mq yr"' , so the time to double the stel- 
lar mass is a significant fraction of the dynamical time of the 
parent GMC, and the main sequence lifetime of a massive 
star. As long as the mass in stars M* < ecMc the cluster does 
not strongly affect the GMC. However, once the stellar mass 
reaches this upper bound, the cluster begins to disrupt the 
cloud, eventually cutting off its own fuel suppl y. The disrup- 
tion t akes less than the GMC dynamical time (IMurrav et alj 
l20IOh . so that further star formation will only increase the 
cluster luminosity by a factor of a few at most. The radiation 
force exerted on the surrounding gas can therefore only be a 
few times larger than the gravitational force. As a result, the 
ejection velocity with which radiation pressure can expel gas 
is of order the escape velocity of the system: 



V2GMci/rc: 



2.2. Size and mass distributions 



(3) 



Massive star clusters are compact. For example. Milky Way 
globular clusters have half light radii r/, ~ 2 - 3 pc, indepen- 
dent of cluster mass, for clusters with masses between lO"* 
and 1O''M0 (see figure[T]i. Observations of young clusters in 
other galaxies al so find rj, ^ 2 p c independent of luminosity or 
Mci (Scheepma ker et al.ll2007h . Younger clusters of t he same 
mass tend to have even sm aller half-light radii, e.g. iLarsenl 
( 12004 and i McCrady & Graham (2007) . This is expected if 
the gas sloughed off by stars as they evolve is also ejected 
from the cluster; the lack of intercluster gas strongly implies 
that the latter is the case. 

Clusters with M^ > 1O^M0 do show an increase in radius 

with increasing mass, of the form r^ ^ mJ (iHasegan et alJ 
2005; Murray 2009); see Figure[T] Numerically, the velocity 
of a cluster is given by 



100 



Md 



lO^Mp 







km s 



(4) 




Fig. 1 . — The present day half hght radii for star clusters. Filled squares 
depict Milky Way globular clusters, filled triangles M 31 globulars, open tri- 
angles correspond to M82 superclusters, filled pentago ns are globular cluster s 
from Cen A. Filled hexagons are ultra-compact dwarfs iHasegan et al.l 120051) : 
IHilkeretalJf2D07l) 

where a = 1/2 for M^ < lO^M© and 0=1/5 for more mas- 
sive systems. As discussed below, this quantity is crucial in 
determining whether a newly formed cluster can expel its sur- 
rounding material above a galactic disk and beyond. 



The star cluster mass function is known to be shallow: 



dN 
dM7, 



ozM-P 



(5) 



where A^ is the number of clusters of mass M^, with /3 « 1 .8. 
This suggests that a few to a few dozen clusters dominate the 
luminosity of a star forming galaxy, and hence certain aspects 
of the ISM dynamics. Direct counts of the number of clus- 
ters that provide half the star formatio n of the Milky Way 
jMurrav & RahmanI 120101) and M82 dMcCradv & GrahamI 
l2007h (see their figure 8) are consistent with this conclusion. 

3. DRIVING WINDS FROM MASSIVE STAR CLUSTERS 

3.1. The disruption of giant molecular clouds 

At early times, i.e. less than a few Myr after a massive 
star cluster forms, the forces acting on the surrounding mate- 
rial include only radiation pressure and shocked stellar winds. 
The first supernovae explode after about 4 Myr. Observa- 
tions show that massive star clusters disrupt GMCs both in the 
Milky Way ( Murray & Rahman 20101) and in nea rby galaxies, 
includi ng the large Magellanic C loud (IOevlll99^ and the An- 
tennae ([Gilbert & Grahamir2007h . As shown by Murray et"aI1 
(l2010h . this disruption is primarily due to the effect of radi- 
ation pressure. While it has been thought that the dynamics 
of gas around star clusters was dominated by shocked stel- 
lar wind pressure, observations of HIl regions in the Milky 
Way have shown that the pressure in such hot gas is equal 
to that of the associated HII (lO'^K) gas (iMcKee et allll984l: 
iHarper-Clark & Murray|2009t iMurrav et al.ll2010l) . and there- 
fore dynamically irrelevant. 



Why don't shocked stellar winds affect the bubble dynam- 
ics? Examination of actively star forming regions in the 
Milky Way shows that the bubble walls are far from uni- 
form; multiple shells, clumps, and 'elephant trunks' (pil- 
lars pointing toward the center of the bubble) are common 
(iMurray & Rahman 2010) . This is not unexpected, given the 
turbulent nature of the ISM; the gas near the cluster is clumpy, 
even before any star formation begins. The fact that the bub- 
ble walls do not actually isolate the bubble interior from the 
region outside the walls explains both the low luminosity of 
diffuse x-ray emission inside the bubbles, and the dynami - 
cal irrelevance of the hot gas (IHarper-Clark & Mur ray 2009). 
However, it is worth stressing that there are clear shell struc- 
tures surroun ding the bubbles in a l l of th e star forming regions 
examined by iRahman & Murravl ( 120101) . and that the bubble 
walls are moving radially outward at 10 -20 km s~'. 

As outlined above, once the radiation pressure of a newly 
formed star cluster balances the gravitational force acting on 
the surrounding gas, the star formation rate slows. The lumi- 
nosity gradually increases, driving gas out of the cluster. As 
the system evolves, a bubble in the interstellar medium forms 
around the cluster. 

For our purpose, we can ignore the fact that the shells are 
clumpy. Since radiation pressure supplies the bulk of the mo- 
mentum absorbed by the gas, the clumpy nature of the ISM 
will not prevent the formation of a bubble: any gas in the 
vicinity will absorb momentum from the radiation field. Some 
very dense clumps will not be appreciably accelerated, but 
the bulk of the gas in a GMC is in a relatively diffuse state 
(n ^ lOOcm"^) and will be pushed outward, piling up in a 
shell, or in multiple shells. 

Over time the bubble expands, so the surrounding shell 
sweeps up more gas, and the shell mass Msh(r) increases with 
increasing radius. The evolution of the optical depth depends 
on the surface density through the shell, Ej/, = Msh(?")/47rr^; 
for a Larson-like GMC density profile p{r) ex r"', Si/,(r) is 
constant. When the shell, or some part of it, breaks out of 
the GMC, near the surface of the gas disk, the growth rate of 
Msh(/") slows and eventually halts. After that, the column of 
the shell or its fragments will decrease with increasing radius. 



iMurray et al.l (1201 Oh studied the r < H behavior of such a 
shell in a GMC centered at the disk midplane. Their model 
included forces due to protostellar jets, shocked stellar winds, 
HII gas pressure, turbulent pressure, gravity and radiation 
pressure. The jets exert an outward force, but are relevant 
only on small scales ~ 1 pc; the shocked stellar winds may 
also be confined, briefly, on these scales; we will refer to the 
combination of the two as the small-scale force, Fs,„. These 
authors found that, on larger scales, the radiation pressure and 
the gravity dominate the dynamics when the star cluster has 

Since the largest GMCs have radii only a factor of two or 
three smaller than the disk scale height, bubbles capable of 
disrupting their parent GMCs will often break out of the disk. 
During the expansion, some shell fragments will be acceler- 
ated in the plane of the gas disk, while others will be accel- 
erated vertically, out of the plane of the disk. In this paper 
we are interested in the latter, which we will assume cover a 
fraction Co of the sky as seen from the center of the cluster 
Before any supernovae erupt in the driving cluster and before 
the shell fragments rise above the disk and become exposed 



to radiation or any supernova-driven hot outflows from other 
stars or clusters in the disk, their equation of motion is given 
by 
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df 
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(6) 



where L^ is the luminosity of the cluster. We argue that the 
breakout of shell fragments from the gas disk is the crucial 
step in driving a superwind. We refer to shell fragments that 
manage to rise above the gas disk as clouds. As we will show 
below, radiative driving, acting alone, can drive clouds above 
the disk scale height in a few Myr, i.e., before the explosion 
of any supernovae from the parent cluster. 

3.2. Scaling relations for ejecting gas clouds above the 
galactic disk 

Here we provide an estimate of the critical gas and star 
formation rate surface densities required to launch a galac- 
tic wind, starting from the assumption that GMC shell frag- 
ment breakout is the limiting step. To do so we first estimate 
the mass and radius of the largest giant molecular clouds in 
a galaxy with a specified circular velocity v^, disk radius R^ 
and gas surface density Sg. The surface density Ec of the 
GMC then allows us to estimate the stellar mass of the largest 
star cluster in the GMC. Given the mass-radius relation for 
star clusters described in section |2] we can then compare the 
ejection velocity of the cluster to the galaxy circular velocity 

Vc- 

We assume that galactic disks initially fragment on the disk 
scale height H ~ {vT/vc)Rd (given by hydrostatic equilib- 
rium), where vj is the larger of the turbulent velocity or the 
sound speed of the gas. The fragments will form gravitation- 
ally bound structures with a mass given by the Toomre mass. 



Mg- 



: TlH^T., 



(7) 



The value of vj required to estimate H can be obtained by 
assuming that the disk is marginally gravitationally stable, so 
that Toomre's Q is of order unity; 



Q- 



TTGRdT. 



1. 



(8) 



g 



It follows that H = ■KGQiRd/vcYY.g. The mass of a large GMC 
is then 

Mo = n'G'Q^(^^\l (9) 

The star formation efficiency in a GMC, i.e. the limit beyond 
which the GMC is disrupted by radiation pressure, (see Eq.|2]i, 
can be expressed as 



ttGc 



ttGc 



eg^ 



2(L/M*) ^ 2(L/M,) 






(10) 



where = H /re sa 2- 3, and the light-to-mass ratio (L/M^.) is 
3000cm- s""* for a lChabrierl (1200 5) initial mass function. The 
characteristic cluster mass is then 



TT^G^C . (Rd 



(11) 



The cluster ejection velocity, as given by eq.[3] reads 



^'ei — ^esc ~ 



-K^G^c^l- Q(j) fRdY. 



(2L/M,)i/2 ^^j V V, 



(12) 



If this velocity is comparable to or larger than v^, then a galac- 
tic scale wind will result. We define the velocity ratio 



C^^ 



2v2 



V, " (2L/M,)i/2 77^ —^ 



(13) 



For C -^ 1 we expect the largest star clusters in a galaxy to 
launch galactic scale winds. For a specified disk radius and 
circular velocity, this defines a critical gas surface density in 
order to launch a galactic wind. 

3.2.1. Critical star formation rate 



We now use ( fT3T l, in conjunction with the lKennicuttI (119981) 
star formation law 



S* = ri 



R. 



(14) 



with rj = 0.017, to find the critical surface density of star for- 
mation required to launch a supergalactic wind. We find 



(2L/M,)i/4 jyC rXvc 



1/4,5/2 



7rGci/4 
Scaling to an L* galaxy. 



(15) 



0.06 



Vc 



200 km s" 



'" ^5kpc\^ , 1 

-^j Moyr-'kpc-'. 



(16) 
From observations of galaxies with and without superwinds, 

iHeckmanI ( 12002) gives t'f" w O.lMgyr-'kpc"". Finally, we 
can express our radiatively-motivated threshold in terms of 
critical star formation rate 
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(17) 



4. MODELING GALACTIC WINDS 



We now extend the model used by iMurrav et al.l ( 120101) to 
describe not only the destruction of a cluster's natal GMC but 
also the late-time evolution of a shell fragment that emerges 
from the gas disk. A number of new effects must be included: 
expansion and varying optical depth of the cloud, the radiation 
pressure contribution from neighboring clusters, ram pressure 
from supernovae, and the underlying dark matter potential. 
We present a dynamical model including these ingredients 
and show that radiation pressure from a massive star cluster 
expels cold gas above the host galaxy disk, where a combina- 
tion of radiation pressure and ram pressure from a hot super- 
nova driven wind then pushes the clouds out to scales reaching 
several tens of kiloparsecs. 

Below we study the dynamics and properties of the gas as 
it travels above the disk. Its equation of motion is described 
by 



(18) 



where the radiation and ram pressure forces Fradir,!) and 
Frad(r,t) include contributions from the parent cluster as well 
as other UV sources distributed over the disk and Fgrav de- 
scribes the gravitational effect of the galaxy and its surround- 
ing dark matter distribution. For simplicity we do not include 
the interactions any hot corona that might occ upy t he galaxy 
halo. We discuss its potential impact in section l672l 

4.1. Launching the gas through a galactic halo 

We have shown that radiation from a single massive star 
cluster can eject clouds from the disk of a galaxy. However, 
neither radiation pressure nor ram pressure from supernova- 
driven hot winds arising from a single isolated cluster can 
drive the clouds to tens of kiloparsecs; doing so requires the 
collective effect of all the clusters in the galaxy, as we now 
show. 

4.1.1. Radiative force 

Once the shell emerges from the galactic disk, it is subject 
to radiation pressure from stars in the disk. The observed sur- 
face brightness of star forming disks follows distributions that 
range from nearly constant forR < R^, where R is the distance 
from the galactic center, to exponential, i.e., Y^uv w exp"*/^'' 
dMartin & Kennicutt 2001; Azzollini et al. 2009). For sim- 
plicity we will work with a constant surface brightness out 
to Rd. If the height above the disk (which we approximate 
as r-H, where r is the distance from the center of the star 
cluster) is smaller than Rj, the disk flux seen by the cloud, 
which quickly dominates that provided by the natal cluster, is 
roughly constant. The radiation pressure then becomes 



if r>RD 



'^rad.cl "•" ^ ^2 



(19) 



where r is the distance from the center of the cluster. At early 
times, the optically thick regime applies and the above rela- 
tion does not depend on k. Once the cloud reaches r > R^, the 
flux seen by the cloud drops as l/r^. 

The radiation pressure force also depends on the size of the 
cloud, which expands perpendicular to the direction of motion 
due to its finite temperature. The rate of expansion perpendic- 
ular to the radius vector will depend on a number of difficult to 
quantify factors, including the pressure of any hot component 
in the disk proper. The clouds are overpressured compared the 
ISM as a whole, since they are subject to the intense radiation 
from the central star cluster. If this pressure is large compared 
to the pressure of the hot gas component, which is the case 
while the cloud is in the gas disk, and possibly for heights of 
order the disk radius, the cloud will effectively expand into a 
vacuum; in this case the size of the cloud in the direction per- 
pendicular to its motion will be given by lj_ = l±o + Cs{t-to), 
where l± is the radius of the cloud in the direction perpendic- 
ular to the radius vector, and c, refers to the sound speed of 
the cloud. 

If, on the other hand, the pressure of the hot gas in the disk 
is large, the cloud volume will vary with height above the disk, 
tracking the pressure in the hot gas. The latter will be in rough 
hydrostatic equilibrium, with a scale height much larger than 
that of the molecular gas disk, so the size of the cloud will 
not vary until it reaches a height of order the scale height of 



the hot gas. At that point the hot gas density will decrease as 
1/r^, as will the density in the cloud. In most of the numerical 
work we present here we assume the cloud is overpressured 
compared to any hot gas; this is the case in the Milky Way, 
but the situation may be different in starburst galaxies. 

If the cloud is at a distance less than R^ above the disk, 
any expansion in the perpendicular radius will give rise to an 
increase in the amount of radiation impinging on the cloud. 
For r > /?d, the 1 /r^ radiative flux decrease will tend to cancel 
the effect of the increased area of the cloud, and the amount 
of light impinging on the cloud will increase less rapidly, or 
not at all. 

The expansion of the cloud is also responsible for a third 
change in the radiative force; the cloud eventually becomes 
optically thin to the continuum emission from the disk. The 
fraction of the incident radiation that the cloud absorbs then 
drops like 1/r^, tracking the decrease in the column density 
and hence optical depth t. If r > Rj, the decrease in incident 
flux tends to cancel the geometric growth of the cloud, and 
the radiative driving, which is proportional to r, decreases as 
r increases. 

4.1.2. Ram pressure from supernovae 

The hot gas from isolated supernovae is likely to be trapped 
in the disk, so that it does not affect clouds above the disk. 
However, supernovae in large clusters will have a different 
fate. As we have just demonstrated, the ISM above large clus- 
ters will be expelled from the disk, opening up a pathway for 
the hot supernova gas to escape. We focus on the large scale 
behaviour of this hot gas, since the bulk of the cold gas is ex- 
pelled from the vicinity of the cluster long before any stars in 
the cluster explode. We assume that the starburst lasts much 
longer then the lifetime of a very massive star (which we take 
to be 4Myr), so that we can average over many cluster life- 
times and use a mean supernova rate. 

Supernovae drive a hot wind. The ram pressure derived 
force on a cloud of cross-section ttIj^ with a drag coefficient 
Cz5 is 

^ram = CDTTl]_ph{Vh - \'cf . (20) 



The kinetic energy of the hot wind is given by 

-MHvl = eLsN. (21) 

where e is the fraction of the supernova luminosity that is 
thermalized to produce a hot phase. It is gene r ally found 
to range from 0.01 to 2 jTheis eFaTI 119921; ICole et alj 
11994'; Padoan et al."1997; 'Th ornton et al.l lI998l). However, 
^Strickland & Heckman (2009) find e > 0.3, possibly reaching 
e = 1, in the local starburst M82. 



Following lStrickland & HeckmanI (120091) . we introduce the 
mass loading parameter /? 



Mh = m 



SN+SW 



:0.2/3M*, 



(22) 



where, by definition, /3 > 1. The factor 0.2 in the second 
equality corresponds to a IChabrierl (l2005h initial mass func- 
tion. For very high values of the mass loading /?, the hot gas 
becomes radiative, cooling on a dyn amical time dSilich et al.l 
120041: [Strickland & HeckmanI l2009b . This effectively Umits 
the range of mass loading, I < /3 < 17. 
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Fig. 2. — The velocity (upper panel) and forces in a model of a IO'Mq cluster embedded in a marginally stable (Q = 1) starforming disk, plotted as a function 
of distance from the cluster center. The time since cluster formation is marked on the velocity curve in the upper panel. The forces shown in the lower panel are; 
the force of gravity (red dashed line) due to the star cluster, the self gravity of the gas, and the gravity of the galactic halo, the force due to protostellar jets and 
possibly shocked stellar winds (cyan dashed line), which is confined to small scales, the radiation pressure force (solid blue Hne) and the ram pressure force from 
a supernova-driven hot outflow (sold green line). The radiation pressure force rises as the cloud emerges from the disk and is exposed to radiation from other 
cluster and stars in the disk; similarly the ram pressure force rises as the cloud enters the large scale hot outflow. In both panels, the vertical dotted lines mark the 
radius of the parent GMC, the scale height of the galactic disk, and the e-folding size of the galactic stellar disk (from left to right). 



Using the above relations and considering that v^- <C v/,, the 
ram pressure force can be expressed as 
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0.2CdM^ X 



The radiation pressure on the same cloud is 

a^ L 

Frad = 7^ X min(ruv, 1) - 
47?- c 



The ratio of the two forces is 
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For typical values: Co = 0.5, e « 0.3, and /3 = 4 and tuv = 1 in 
the optically thick regime, this ratio is '^ 1 . At large distances, 
of order 10 kpc, when the clouds are optically thin, the ram 
pressure force will start to dominate. 

4.1.3. Gravitational force 

While in the disk, the shell takes part in the rotational mo- 
tion of the galaxy, so that the force exerted by the galaxy on 
the GMC is cancelled by the centrifugal force due to rotation. 
However, as it rises above the disk, the cloud will lose the 
centrifugal support it enjoyed while pursuing its circular path 
around the galaxy. We take this into account by using 
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if r < Rd- This assumes that the GMC lies at a distance ~ R^ 
from the galactic center 

We model the overall gravitational effect of the galaxy and 
its dark matter distribution by using a singular isothermal 
sphere model for the galactic halo: 



'igrav ■ 



(27) 



This accurately reflects observed rotation curves of spiral 
galaxies out to r w 50kpc, e.g., iCasertano & van GorkomI 
(|19 91). It probably overestimates the force of gravity acting 
on the cloud at large radii (approaching the virial radius of the 
galaxy), where lensing measurements suggest that an NFW 
profile is a better fit, e.g.. lMandelbaum et al.l ( 120061) . 

4.2. Numerical results 

Given the forces introduced above, we now solve the 1-D 
equation of motion and present results aimed at characterizing 
the fate of the shells described above. 

In Figure |2] we consider the case of an M82-like galaxy, 
i.e. with a circular velocity v^ = llOkms"' and a luminos- 
ity L = 5 x 1O"^L0 . For such a galaxy, our idealized model 
(Eq.|9]l predicts a characteristic c luster mass of about 10^M(7) , 
in agreement with observations dMcCradv & Grahamll2007l) . 
As indicated in section 13.21 the ( parameter of such a sys- 
tem is greater than one, indicating that the star cluster can 



drive a shell of gas above the plane. Our numerical calcula- 
tion confirms the analytic result. We now describe the detailed 
behavior of the shell as a function of time and radius. The bot- 
tom panel of the Figure presents the amplitude of the relevant 
forces as a function of radius: 



• small scale force (F,,,,, cyan dashed line): it represents 
the momentum per unit time deposited by proto-stellar 
jets, outflows and possibly shocked stellar winds on 
s mall scales. The de tails of this force are described 
in lMurray et al.l (120101) . Numerous observations of ex- 
pandi ng bubbles in the Milky W av (iMurrav & Rahman 
20101) and in n earbv galaxies (IWhitmore et al.l U 
Gilbert & Grah am 2007) demonstrate that such 
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out- 
flows start before any supernovae occur The amplitude 
of this force quickly decreases, on a timescale of about 
0.1 Myr. 

• Gravity (red dashed line): At early times the gravity 
is dominated by that exerted by the star cluster on the 
nascent shell; for r < Rqmc this force drops like l/r-. 
At larger radii, H < r < Rj, the mass in the shell is 
constant, but the cloud is losing the rotational support it 
enjoyed in the disk, so the effective gravitational force 
increases. At still larger radii the force is dominated 
by the interaction between the cloud and the combined 
gravity of the stars and dark matter in the galaxy proper, 
sothatg(r) = \'^/r. 

• Radiation pressure (blue solid line): initially, the radia- 
tion pressure force originates only from the parent star 
cluster, and the shell is optically thick to far infrared 
photons. As a result, this force initially scales as l/n, 
but the expansion of the cloud gradually reduces the 
far IR optical depth below unity, and the net radiation 
pressure force becomes roughly constant. The radiation 
pressure from the parent cluster lasts only for a period 
of about 8 Myr, i.e. the time for the cluster luminos- 
ity to decrease by a factor of two. If the shell does not 
reach an altitude greater than the height of the disk dur- 
ing that time, it will fall back toward the mid plane of 
the galaxy. If the cloud succeeds in rising above the 
disk, it then feels the radiation pressure from the neigh- 
boring clusters, which increases until r = Rd. Above this 
height, the radiative flux decreases like 1/r^. The frac- 
tion of the luminosity impinging on the cloud is given 
by the area of the cloud divided by r^. On those scales, 
the size of the cloud increases roughly linearly with 
time, but the radius increases more rapidly. As a result, 
the momentum deposition decreases. Finally, due to its 
continuing expansion, the cloud becomes optically thin 
at about half a kpc, so Fmd drops like l/p- beyond this 
point. 

• Ram pressure (green solid line): we have already noted 
that observations of the Milky Way and the LMC show 
that stellar winds are not dynamically important, so the 
ram pressure force is gradually ramped up when the 
cool gas emerges from the disk, around 60pc in figure 
121 It is comparable to the radiation pressure, with the 
exact ratio of the two forces varying with e and /?, as in 
eqn. ( |25] |. The ratio changes when the cloud becomes 
optically thin to UV radiation, around r ~ 0.5 kpc. At 



1000 



100 



5x10' 


] 1 


1x10' 




3x10' 




1x10'^ 




3x10= 




1x10= 


_/ 



0.001 




100 



Fig . 3 . — Outflows launched by clusters in a M82-like galaxy. The different 
colors denote clusters with different masses, ranging from lO' to 3 X IO'Mq. 
We note that clusters with M < 10* Mq only trigger galactic fountains. 

large radii, the ram pressure is the dominant force, caus- 
ing the cloud to accelerate slowly beyond r ~ 50kpc. 
Because the area of the cloud continues to increase with 
increasing r, but not as rapidly as r^, the ram pressure 
force decreases with increasing r. 

From an observational point of view, we note the following 
results: clouds escape the disk in less than 3 Myr, before any 
supernovae explode in the cluster They reach velocities of 
several hundred kilometers per second at a radius of 10 kpc, 
after a time ^ 10 Myr. The clouds reach a distance of 50 kpc 
after- 100 Myr. 

4.2.1. Star-forming galaxy: varying cluster mass 

As shown in Eq. [13] for a given galaxy the main parameter 
defining the fate of gas clouds ejected by radiation pressure is 
the cluster mass — only sufficiently massive clusters can trans- 
fer enough momentum to eject a shell of gas. To demonstrate 
this result numerically, we use the M82-like galaxy parame- 
ters described above, and solve the equation of motion for a 
range of cluster masses, from 10^ Mq up to IO^Mq. The re- 
sults are presented in Figure [3] We can see that, in a halo 
of mass Ml, « lO^'M©, clusters with M = 10^ Mq will launch 
galactic fountains rather than galactic winds. More massive 
clusters launch winds to larger distances from the center of 
the galaxy, with clusters above M = 10^ Mq expelling gas from 
the halo. 

4.2.2. Quiescent galaxies: varying cluster mass 

We now consider the case of a quiescent galaxy: The Milky 
Way. As done above, we study the evolution of the material 
surrounding young clusters. In such a system, v^ = 220km s~' , 
Alg = 2 X 10^ M0 and Rd = ^ kpc. The results are shown in 
Figure|4] As a result of the Toomre criterion (see Eq.|9]), typi- 
cal star clusters are expected to have M^i =10"*- 10''Mq. The 
trajectories of shell fragments expelled by such star clusters 
are shown in solid fine. Those of hypothetically more massive 
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Fig. 4. — Outflows launched by clusters in a Milky Way like galaxy. The 
dashed lines represent the trajectories of material launched by clusters with 
a mass exceeding that allowed by the Toomre criterion. Our model shows 
that such a galaxy can only trigg er g alactic fountains, in agreement with the 
scaling relation introduced in Eg. 1131 

clusters are represented with dashed lines. As can be seen, a 
galaxy like the Milky Way is able to launch galactic foun- 
tains to heights of order a kiloparsec. However, even clusters 
with M» = IO^Mq, similar to the most massive clusters in the 
Milky Way, do not drive large scale outflows. Such a behavior 
is in agreement with the scaling relation presented in Eq. [T3] 
It reflects the fact that in such a galaxy, the initial kick result- 
ing from radiation pressure is not strong enough to trigger a 
large-scale outflow. 

4.2.3. Starburst galaxies: varying cloud mass 

We now model the evolution of shells of various masses in 
a typical z = 2 starburst galaxy. We parametrize such a system 
withvc=100kms"', Mg = 2x IO'^Mq, 7?^ = 5 kpc andL = 
4 X 10" L0. We consider the case of a Lyman break galaxy 
for which most of the luminosity comes out in the optical. 
Given Eq. |9] the star cluster mass expected for such a galaxy 
is about 3 x 10^ Mq. This stellar mass will be distributed in 
a number of clumps, reducing the self gravity. To take this 
effect into account, we increase the timescale over which the 
radiation pressure from the parent cluster is on by a factor 
four 

Here, as an illustration, we use do not include the ram pres- 
sure contribution from supernovae and consider only radia- 
tion pressure . We find that, for such a galaxy, radiation 
pressure can expel the entire shell with a typical velocity 
V ~ 1 00 km s~' . We also consider the case of lower mass shell 
(or shell fragments): these can easily reach velocities of order 
1 000km s~' and might be related to the tails of blue-shifted 
self absorption observed in star forming galaxies (Weiner et 
al. 2009, Steidel et al. 2010). 

5. CLOUD AND GLOBAL WIND PROPERTIES 

We have shown that radiation pressure from a single mas- 
sive cluster can drive cool gas above the disk of a galaxy. Sub- 
sequently, the radiation and supernova supplied ram pressure 



Fig. 5. — Outflows launched by a z ~ 2 starburst galaxy. Small shell frag- 
ments can reached outflow velocities of order 1 000 km s"' . 



from the collection of star clusters in the disk, can drive the 
cool gas to radii of tens or even hundreds of kiloparsecs. This 
outflow will be in photoionization equilibrium, with a temper- 
ature of order 10,000K and an ionization parameter of order 
U ^ lO"'*. It will be dusty, and it will have the bulk of its 
ions in a moderately low ionization state. The columns seen 
in absorption will range fiom Nh ~ 10^^ cm"^ downward, de- 
creasing with increasing radius and velocity. The cover factor 
Cf of the individual outflows will also increase with increas- 
ing radius and velocity. 

The column densities of the clouds that are launched are 
initially those of the parent GMCs, i.e., Nh ^ 10^^ cm"^ and 
upwards. The column decreases roughly as ScCrc/r)^ once 
the shell has disrupted the parent GMC, an approximate rela- 
tion that holds exactly if the radial velocity is constant at large 
radii. 

The bulk of the gas in starburst galaxies is in molecular 
clouds (unlike quiescent galaxies such as the Milky Way, 
which has a substantial mass of atomic gas). It follows that 
most of the mass in the bubble shell (and later in the shell 
fragments in the wind) originates in the natal GMC. However, 
it is probably true that some gas lies either above or below the 
GMCs in which star formation occurs. This overlying gas 
will slow or possibly halt the bubble shell, if it has a column 
similar to that of the GMC. However, we note that star clus- 
ters in the Milky Way are seen to have a dispersion in their 
distance from the dis k midplane that is compar able to the gas 
scale height H, e.g., (iMurrav & Rahmanll2OI0l) . As a result it 
seems likely that a substantial fraction of clusters above the 
critical cluster mass will drive winds out of one side of the 
disk. 

It seems likely to us that clusters in particularly large 
GMCs, with GMC radii comparable to H such as those seen 
in chain or clump galaxies, can expel winds from both sides of 
the galactic disk. Our one dimensional models cannot address 
this question without further information on the distribution of 
gas outside GMCs. More detailed modelling will be needed 
to decide this question. 



5.1. Cloud properties 

If a wind is launched from a cluster, the column density 
through the wind will be given roughly by 



NH{r)^Ng 



H 



(28) 



where Ng = Yig/nip and /3 « 1-2. The power law index /? 
depends on the details of the model, in particular on the initial 
size of the cloud, how v(r) varies, and on whether the cloud 
is pressure confined by hot gas. The initial column density is 
simply that of the disk (when the GMC has expanded to the 
size of the disk scale height). As already noted, the gas will be 
in photoionization equilibrium, and illuminated by thousands 
of O stars, so its temperature will be of order 10"* K. Since the 
shell or cloud is being accelerated, the width of the shell in 
the direction of the acceleration will be 
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where the acceleration g w v-/r. Since v w v^ ~ v^, we have 
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Combining equations (l28T l and (|30] |, the number density is 
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where we have taken /3 = 2 and v^ = 100 km s~'. For a star 
formation rate of SM© yr"', this corresponds to an ionization 
parameter 
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where Q is the number of ionizing photons emitted by the 
galaxy per second; we have scaled to a value appropriate to 
the assumed star formation rate. If /3 < 2, the r dependence 
is weaker, e.g. /3 = 1 implies U = const. We hav e run the 
photoionization code CLOUDY jFerland et aljl998l) to deter- 
mine the ionization state of the gas. The bulk of magnesium 
is in the form of Mgll for U = IQT^, while a fraction of or- 
der 0.1 of sodium is neutral. At very large r the background 
UV flux will dominate that from the galaxy, and the ionization 
fractions of both ions will drop, even for the case /3 = 1 . 

The column densities, densities, sizes, and ionization pa- 
rameters just given are those of individual cloud in the outflow 
from a single cluster Observations of expanding gas around 
clusters in the Milky Way make it clear that while shell-like 
structures do form, they are not very regular, nor is all the 
material at the same radius from the cluster. Physically, it 
is clear that higher column density gas will accelerate less 
rapidly than lower column density gas, leading to clouds at a 
range of radii. Further, it is likely that material at the fringes 
of the clouds will be ablated, as is seen in star forming regions 
in the Milky Way. When seen in the large, the outflows from 



a cluster will appear quasi-continuous; if the flows from mul- 
tiple clusters cannot be resolved, the outflow from the galaxy 
will appear even more continuous. 

5.2. Absorption covering factor 

We note that the outflow will not be seen as such until the 
continuum (dust) optical depth drops below unity; only then 
may absorption lines be seen against the cluster or galactic 
light. The corresponding column is Nh ~ 10^' cm"^, depend- 
ing on the wavelength of the absorption line used to trace the 
outflow. 

While the continuum optical depth drops below unity at 
fairly small radii, the Mgll resonance transition at 2900 A is 
optically thick out to radii as large as 100 kpc. The dynamical 
time to lOOkpc is ^ lOOMyr. Hence we expect blue shifted 
absorption to be seen in post-starburst galaxies as well as in 
active starbursts. 

Individual GMCs and their cluster progeny occupy a small 
fraction of the disk area, so individual outflows have a small 
disk covering factor 



C 



fd- 






(33) 



This covering factor is an increasing function of cloud veloc- 
ity V for small distances r, since initially both v{r) and a{r) are 
increasing functions of r. When the clumps enter the large 
scale hot outflow, their radii may decrease as they come into 
pressure equilibrium, but after that they will resume their ex- 
pansion. However, the summed cover factor of all the clouds 
will be a significant fraction of unity. 

The ratio of disk scale height to disk radius is related to the 
star formation rate of the galaxy. For high star formation rates, 
it can reach (H/Ri)^ w 0.1-0.3, but is smaller for most star- 
burst galaxies. However, clusters above and below the disk 
midplane can drive gas from either side of the disk (and pos- 
sibly from both sides, as just noted), so that sightlines from 
most directions will encounter the outflow from one or more 
clusters. 

5.3. Mass loss rates 



iMurray et"an (120 1C<) show that the fraction of gas in a GMC 
that ends up in stars is given roughly by 
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interpolating between their optically thin and optically thick 
expressions. The balance of the GMC is either returned to the 
interstellar medium or driven out of the galaxy as a wind. We 
have argued i n this paper that for systems above the thresh- 
old given in ^3.2.11 one quarter to one half of the GMC is 
launched as a wind. The ratio of mass launched out of the 
disk to mass locked up in stars is very roughly 



M„ _ 0.25gcm-2 + S]G 



(35) 



This predicts a very large ratio of mass lost in a wind com- 
pared to stellar mass in low Sc systems. Since we require 
that the galaxy be above the threshold star formation rate, this 
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high mass loss ratio applies to galaxies with low circular ve- 
locities. 

In addition to this local criterion, there is a global limit to 
the mass loss rate driven by radiation pressure, namely that 
the rate of momentum carried out by the wind is no larger 
than the rate of momentum supplied by starlight. This limit 

is 
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The maximum wind mass loss rate per star formation rate for 
such a radiatively driven outflow is 
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These expressions assume that the bulk of the star formation 
occurs in clusters above the critical mass necessary to launch 
galactic winds. As noted in section IZ21 the mass function of 
star clusters is shallow so that this criterion will be satisfied 
for galaxies with sufficiently high star formation rates. 

6. DISCUSSION 

6.1. Comparison to other launching mechanisms in star 
clusters 

6.1.1. Ram pressure from supernovae and cloud survival 

There is little doubt that hot gas from supernovae pushes 
cool gas out of galaxies. However, it has been argued here 
and elsewhere that supernovae do not push cool gas out of 
GMCs. There is good evidence in the Milky Way that GMCs 
are disrupted well before any supernovae explode in th em 
iHarper-Clark & Murray! (120091) : iMurrav & RahmanI (120101) . 

The crucial point, however, is that cool gas entrained in a 
hot flow on GMC scales will not survive at a low tempera- 
ture long enough to escape the galaxy. As noted in the in- 
troduction, clouds caught up in a hot outflow are compressed 
on a cloud crossing time, t^c = x^^^'^hb- For a cloud with 
T = lOOOK at a distance r ss lOpc from its natal cluster, we 
find 

(39) 
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where we have scaled to the temperature Th = 5 x 10 ^ K of the 
hot flow seen in M82 (IStrickland & Heckmanll2009h . 

Numerous simulations of interactions of hot density gas 
with cool (r < IO^'K) gas demonstrate that the clouds are 
ra pidly destroyed, i .e., i n ^ 3 - 5tcc, in t he sim ulations 
of lKleinetail (11994 and iPoludnenko et al.1 (120021) . These 
early simulations employed an adiabatic equation of state, 
so that the temperature of the cool cloud increased as it 
was compressed. Simulations which include radiative cool- 
ing find that clouds live longer, ~ lOf^e, re aching veloci- 
ties o f several hundred kilometers per second (iMellema et al.l 
I200I ICooper et al.1 l2009t iPittard et al.1 I20T0I) . For exam- 
ple. lCooper et al.l (l2009h emplov high resolution 3D radiative 
models, finding that a substantial fraction of the gas in a cloud 



is lost, via Kelvin-Helmholtz and Rayleigh-Taylor instabili- 
ties, to the hot flow by the time the cloud has traveled 75 pc 
(about a million years), consistent with the simple analytic 
estimate of ^ lOfcc- 

The short Ufe times of such ram pr essure driven clouds 
was stressed by iMarcolini et al.l (l2005h . who noted another 
form of destruction, that of thermal conduction. They find 
that conductive heating suppresses the Kelvin-Helmholtz and 
Rayleigh-Taylor instabilities that disrupt the clouds in non- 
conductive simulations. However, the conduction also causes 
rapid mass loss from the clouds, l eading to life t imes ~ 1 - 
lOMyr. As in the simulations of ICooper et al.l (l2008h . the 
clouds are essentially completely disrupted by the time the 
reach r w 1 kpc, even in the most favorable simulations. 

Observations of the solar wind have shown that the conduc- 
tivity is similar to the classical Spitzer value, e.g.. lSalem et al.l 
( l2003i) , modified to account for saturation at the free- 
streaming heat flux. This, combined with the rough agreement 
between the analytically estimated cloud lifetimes against 
conduction and t hose f ound by the numerical simulations of 
IMarcolini et al.l (120051) . strongly suggest that ram pressure 
driven cold gas clouds have lifetimes of order ~ 1 Myr, simply 
due to heat conduction, independent of the effects of Kelvin- 
Helmholtz instabilities. 

To summarize, the destruction time fdest ^ 20fcc of cold 
clouds entrained in a hot outflow is given by 
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Assuming a cloud velocity of 300km s , the cloud is de- 
stroyed after traveling a distance 
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This is a rather generous estimate, since it assumes the upper 
limit for estimates of the cloud destruction time, and assumes 
that the clouds are accelerated instantaneously. The sizes of 
starburst disks range from R^Ri 0.3 -I kpc or larger, implying 
that the hot gas density does not begin to drop until r = 0.3- 
Ikpc. 

We conclude that cold clouds entrained in hot flows on 
GMC scales do not survive beyond distances of '--^ 1 kpc. 

Radiatively driven cold gas clouds escape their natal clus- 
ters well before any supernovae explode, so they are not sub- 
ject to the ill effects of hot gas, at least until they reach the 
surface of their galactic disk. If the hot gas fills the volume 
above the disk surface, the cold clouds will be subject to both 
destruction by Kelvin-Helmholtz instabilities and conductive 
evaporation, but at a much lower rate than the clouds simu- 
lated above. The slower evaporation results from the fact that 
the conductivity is saturated; in that case the mass loss rate 

A^cond ^ «/, , where «;, is the number density of the hot gas. 
The longer Kelvin-Helmholtz times follow from the larger 
length scale a of the clouds, which grows as the clouds move 
away from the launching cluster 

This leads to a clear observational distinction between the 
two types of driving. Small clouds driven by hot gas have 
lifetimes of order 1 Myr or shorter in most simulations; in the 
most optimistic case the Hfetime may reach ~ 3 Myr. The cold 
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clouds survive to r < 1 kpc from the disk (assuming a rather 
generous instantaneous v = 300km s~ ). 

In contrast, radiatively driven clouds are not subject to such 
rapid initial mass loss. They will travel for ^ 4Myr, reach- 
ing distances 300pc before the hot gas from the first cluster 
supernovae explode, and about twice that before the bulk of 
the supernovae in the cluster explode. The cold clouds will 
be traveling at ~ 100km s~' or more when the hot gas reaches 
them, and they will have sizes of order ^ 50 pc, assuming only 
that they expand at their sound speed (initially the clouds may 
expand more rapidly, as the radiation field diverges). As a re- 
sult, radiatively driven clouds will survive to reach r ^ 50 kpc 
or more. Cool outflows extending to r ^ 1 kpc from the host 
disk are a clear signature of a radiatively driven outflow. 

A second test that can distinguish between radiative driv- 
ing and supernova driving is even more direct. Individual star 
clusters are known to drive winds. If radiative driving is im- 
portant, winds will be seen emanating from clusters younger 
than 4Myr. Since no stars have lifetimes shorter than this, su- 
pernova driving cannot account for any winds seen emerging 
from such young clusters. 

6.1.2. Shocked stellar winds 

Other mechanisms have been proposed to drive bubbles in 
the Milky Way and winds in other nearby galaxies. For ex- 
ample, shocked stellar winds from O stars produce lO^K to 
10^ K gas seen via X-ray emission. This gas was suggested to 
be the driving force behind the bubbles observed in the Milky 
Way by lCastor et al.l (Il975h . However, on closer examination, 
these models fail to reproduce the observations. T he predicte d 
bubble radii are much larger than observed, e.g. (IOeylll99a) . 
as is the diffuse X-ray emission from the shocked gas. Thus 
the shocked stellar winds either cool in place without emitting 
X-rays (McKee et al. 1984) or escape from the bubble thro ugh 
holes in the bubble wall t H arper-Clark & Murray! 120091) . In 
either case they do not affect the dynamics of the bubble. 

6.1.3. Other hybrid models 



Both lMurray et all (120051) and lNafli & STIS (120091) discuss a 
combination of radiative and ram pressure driving for galac- 
tic superwinds, but neither paper addresses the fact that most 
stars form in clusters, nor do they mention the short lifetimes 
of ram pressure driven outflows originating at small radii. The 
latter implies that cold outflows are restricted to a few kilopar- 
secs around the host galaxy, in contrast to the prediction made 
here, that cold gas survives to 10- 100 kpc. 

We find that clusters are crucial for launching winds; the 
radiation from massive clusters blows holes in the local ISM, 
starting the cold gas on its way, and allowing for the hot gas 
from subsequent supernovae to escape readily. This hot gas 
can then drive cold gas, both that from its own natal cluster, 
and cold gas from other, older clusters to large radii. 

We have stressed that launching cold gas above the disk is 
the crucial step in driving a superwind. We used this point to 
estimate the critical star formation rate surface density. Be- 
cause the radiative force is less destructive than ram pressure 
driving by hot gas, and because the clouds grow in size as 
they move away from their launching points, the clouds can 
survive the any subsequent blast from hot gas. 



6.2. The hot corona 

We have neglected any hot gas which might reside in the 
galaxy halo, and which could potentially halt the outflow of 
cold gas clouds within the halo. This is drag (or destruction) 
is relevant for any proposed launching mechanism. 

If the halo co ntained the cosmic abundance of baryons, 
Mh ~ 0.173M/, dPunklev et al.1 120091) . the mean column of 
hot gas would be Nh ^ 3 x lO'^cm"-. Since this is much 
larger than the cloud column density at the virial radius, 
the hot gas would either stop or shred the outgoing clouds. 
However, the hot gas content of cluster and group size ha- 
los has bee n measured to be smaller than the cosmic abun- 
dance, e.g., iDai et al.l ( 120091) . These authors measure a hot 
gas fraction fg that ranges from fg w 0.08 ± 0.02 for galaxy 
halo masses M,,=6x IO'^^Mq down to fg = 0.013 ±0.005 for 
Ml, = 2.6 X IO^^Mq. Written as a ratio to the cosmic baryon 
abundance, this is fg/fc = 0.46 to 0.075. It is believed that 
the hot gas fract ion in galaxies (M/, < IO^^Mq) is yet smaller 
dBregmanl [20071) . If we take fg/fc = 0.1, flien cold clouds 
reach r ^ 50- 100 kpc before they run into their own column 
of hot halo gas. 

7. CONCLUSIONS 

We have shown that radiation pressure from massive stel- 
lar clusters can drive outflows with velocities v ~ a few x 
100 km s~' to distances ~ 50- 100 kpc from their host galax- 
ies. The outflows consist of cold {^ lO'* K) and dusty 
clouds with characteristic column densities in the range A^^ ~ 
10^' cm"^ near the galactic disk, ranging down to Nh ~ 
10"* cm"^ at r ~ 50kpc. The gas in the clouds is in a low 
ionization state. Once a cloud is above the disk of a galaxy, 
it experiences the radiation pressure from nearby clusters as 
well as ram pressure from supernovae. 

This solves a major problem of hot outflow-driven cold 
clouds, namely their short lifetimes. Since radiatively driven 
clouds need not be in contact with (relatively) dense hot gas 
in order to be driven out of the galaxy, they are not destroyed 
by Kelvin-Helmholtz and/or conduction. 

We have argued that individual clusters can drive outflows 
above the disk. If the disk is sufficiently luminous, the disk 
light can then drive the outflow to distances of tens of kilo- 
parsecs. This evades the need for a global "blowout" of the 
interstellar medium of (the central portions of) the galaxy in 
order to drive a large scale wind. We note that ULIRG are 
clearly not "blown-out" and yet have winds, so this is an at- 
tractive feature of cluster driven winds. 

Since the gas is launched from individual star clusters 
roughly vertically from the disk, it will retain the rotational 
motion of the disk. If the launching cluster is at large disk 
radius, measurements of the rotational velocity will find gas 
at large radii and large azimuthal velocities. The mass loss 
rates from cluster driven winds are comparable to the star for- 
mation rate in massive galaxies, and substantially larger than 
the star formation rate in low mass galaxies. In both types of 
objects, the terminal velocity of the wind is similar to the rota- 
tional velocity of the galaxy. While the winds may or may not 
reach the escape velocity, they can nevertheless reach radii in 
excess of the host galaxy's virial radius. 

Since the cluster escape velocity must be comparable to 
the escape velocity from the galaxy to launch a flow above 
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the disk, and because the maximum cluster mass is related 
to the gas surface density and hence star formation rate, 
there is a critical star formation rate (per unit area) E» « 

— 1 — "^ 

0. 1 Mq yr kpc " to launch a wind radiatively. This is in good 
agreement with the threshold inferred from observations. 

This paper has employed simply analytic and one dimen- 
sional hydrodynamic models of cluster driven winds to ar- 
gue that such winds are crucial to the formation of galactic 



superwinds. The results provide strong motivation to carry 
out more realistic three dimensional radiative hydrodynamic 
modeling. 
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